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Abstract 

The structural relaxations of a dense, binary mixture of charged hard spheres are studied 
using the Mode Coupling Theory (MCT). Qualitative differences to non-ionic systems are 
'\q- shown to result from the long-range Coulomb interaction and charge ordering in dense 

molten salts. The presented non-equilibrium results are determined by the equilibrium 
structure, which is input using the well studied Mean Spherical Approximation. 



Introduction 



The equilibrium structure of ionic liquids is strongly affected by the long-range nature 
of the Coulomb interaction [1]. One aspect is the screening of external charges familiar from 
the Debye-Hiickel theory of ionic solutions. Another aspect are oscillations in the charge 
density around a given ion. These oscillations result from the competition between local 
^j- ■ charge neutrality and the excluded volume restriction due to the finite diameter of the ions. 

Also the dynamics of equilibrium ionic liquids in the hydrodynamic regime differs from the 
one of uncharged mixtures. Coulombic restoring forces lead to a non-diffusive and non- 
propagating relaxation of charge fluctuations [1]. The MCT determines the slow structural 
relaxations of dense (supercooled) liquids from their equilibrium structure [2] . In this paper 
we discuss the most salient features of these results arising from the long-range Coulombic 
interactions. These will be seen to be directly connected to the familiar effects in the static 
. structure mentioned above: screening and charge ordering. 

: 

Theory 

One of the most simple models of an molten salt is a binary mixture of hard spheres 
with radii di and charges z iy i — 1, 2. Global charge neutrality fixes zig± + z 2 g 2 = 0, where 
denotes the density of species %. The mean spherical approximation [3] gives a satisfactory 
description of the equilibrium structure of this system, which depends on the parameters 
d\/d2, Z\j Zii packing fraction ip and coupling constant T. The packing fraction is the ratio of 
volume occupied by spheres to the total volume. T is a generalized inverse Debye screening 
length and a measure of the Coulomb interaction compared to the thermal energy. At the 
values of these parameters chosen in our study and collected in table I one observes a liquid 
to glass transition in the MCT equations. This transition is the topic of the work reported. 
It is worth noting that the density at the transition can be chosen to be lower than in 
the uncharged system [4]. Obviously, the charges increase the interactions of the particles. 
Asymmetric parameters (di/d 2 ^ 1) were first chosen in order to study experimentally more 
realistic non-symmetric salts. The charge asymmetry zi/z 2 , however, was then adjusted to 
obtain the value A = 0.85 for the exponent parameter A; see the discussion. The MCT for 
binary mixtures formulates a closed set of equations for the time and wave vector dependent 
density fluctuation functions, F^ b {t) = jj(dQq*(t)6gq(0)) [5,2]. Of particular interest are the 
fluctuations of the total- or mass-, g n = g 1 + g 2 , and the charge-density, g c = z\g x + z 2 g 2 . 



Table I: Molten salt equilibrium parameters and results from MCT calculation. 
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The initial values, F ' (i = 0) = 5" , are the static structure factors which are the only 
input determining the F°- b (t) via the MCT equations. 

MCT shows that these equations exhibit bifurcations which are identified as idealized 
liquid to glass transitions. Close to these transitions the long time dynamics is predicted to 
follow from 

q 2 F q (t) = S q {M q (t)S q - j t jf* dt'M q (t - t')F q (t')} , (1) 

where matrix notation is used and the memory functions are quadratic polynomials in the 
F q (t) correlators with coefficients determined by the static structure factors S q [2,5]. Ther- 
mally activated transport is neglected in (1) leading to its breakdown at low temperatures. 
Only aspects which are not affected by this simplification will be discussed in this article. 

Results 

In order to screen an external charge the charge structure factor S q c has to vanish for 
small q, S q c oc (gA^) 2 , where is the Debye screening length [1,3]. This leads to a 
decoupling of the MCT equations (1) for small wave vector. Whereas the mass-density 
fluctuations are determined by a frequency dependent longitudinal viscosity N l (z = cu + ie), 

F™(z)/S™^ ~\ for q^0, where N l q (z) = \m™(z) , (2) 

the charge fluctuations couple to the generalized conductivity a(z) 

F^(z)/S^ c -> ~ l . . . for q -> 0, where a(z) = lim - f € . . . (3) 

q w/ q z + Ama(z) H v 1 <?^o 4vr S c q c M^(z) v 1 

These equations simplify further in two frequency windows reached close to a transition. 
From the well known MCT results let us only mention the von Schweidler decay which 
describes the onset of the a-relaxation, i.e. the final decay into equilibrium [2]: 

Ff(t)/Sf = ff - hf (t/r) b for intermediate times . (4) 

The von Schweidler exponent b and the exponent 7 determining the increase of the a- 
relaxation time r are functions of the exponent parameter A and uniquely determined at 
the chosen transition. Their values are included in table I. Eq. (4) shows that the density 
fluctuations exhibit a two-step relaxation. The amplitudes f q b of the final or a-relaxation 
is smaller than unity. The relaxation is non-exponential in general and the ct-relaxation 
times Tg b depend sensitively on temperature or density and on wave vector by a multi- 
plicative factor, which roug hly equals (/f /hf ) 1 / 6 . The mass-density correlation functions 
qualitatively agree with the results obtained for neutral one-component liquids [2,6]. Fig. 
1 shows the a-amplitudes, f q n , which describe the frozen-in mass-density structure at the 
glass transition. 




Figure 1: Mass-density a-amplitude critical amplitude h™ n and structure factor S™. 

The known local packing on shells separated by the mean average interparticle spacing, 
a m g n /4.4, is seen in /™ n as a consequence of the one in S^ n [2,6]. Local neutrality and 
excluded volume effects lead to charge ordering and a prominent peak in the charge structure 
factor, S™. The average spacing between ion-shells of equal sign is larger than the average 
particle spacing resulting in the peak in S^ c to lie at q c , where q c < q n . The charge-density 
fluctuations f^ c , which are arrested at the transition, reflect this ordering [4]; see Fig. 2. 




Figure 2: Charge-density a-amplitude f£ c , critical amplitude and structure factor S\ 
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The different mass- and charge-density ordering also leads to specific variations in the 
wave vector dependent prefactors of the a-relaxation times r" n and r^ c . Fig. 3 shows 
these variations as estimated from the peak positions in the corresponding susceptibilities, 
T~uj max = 1. A g-dependent slowing down at the maximal amplitudes of the a-process is 
observed. Superficially, this mimics the known De Gennes narrowing, as the T q vary in 
phase with the static structure factors [1]. However, this correlation only holds, because 
the a-relaxation amplitudes, f q , vary in phase and the critical amplitudes, h q , vary out of 
phase with the structure factors, S q [2]; see Figs. 1 and 2. 

The results presented thus far specify the structural relaxations on length scales of the 
order of interparticle distances. Analyzing Eq. (1) more closely, it is seen that it also is 
the static structure on these length scales which determines the results. The quantitative 
results therefore depend on the appropriateness of the underlying microscopic model, i.e. 
the liquid of charged hard spheres with the choice of parameters. From Eqs. (1,3) one 
can also obtain macroscopic transport coefficients like the conductivity a and the dielectric 
constant e' of the ionic melt. Their values are included in table I. The frequency dependent 
conductivity determines the dielectric "constant" via 

e(z) = 1 + 47ri^ . (5) 

The liquid molten salt is characterized by a conductivity and dielectric constant measured 
at low frequencies, wr C 1. In the idealized glassy state the particles are arrested and ionic 
transport over macroscopic distances is not possible. The a-relaxation time r diverges and 
the conductivity vanishes. Eq. (5) then results in a dielectric constant which can also 
be observed in the liquid state at high frequencies, out ^> 1. Fig. 4 shows the dispersion 
of the dielectric constant versus ujt. In the same plot, the conductivity crosses over from 
its low frequency value to a power law behavior, a(ur ^> 1) oc (cjt) 1-6 , at intermediate 
frequencies. 
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Figure 3: Wave vector dependent factors of the a-relaxation times r" n /r and r c q c jr. 




Figure 4: Scaled conductivity or and real part of the dielectric constant e versus lot. 



Discussion 

The mixed salt CKN [40% Ca(N0 3 ) 2 - 60% KN0 3 ] is a well studied glassforming melt; 
see [7,8] and references therein. Its dynamics in a wide temperature range has been studied 
by neutron [7] and by depolarized light scattering [8]. In the latter work it was found that 
the exponent parameter A ~ 0.85 describes the dynamics in an intermediate time window. 
These findings lead to the choice of parameters in our model reproducing this A-value. This 
assures that the asymptotic dynamics in the intermediate time window is described by the 
correct asymptotic master function exhibiting, for example, the von Schweidler asymptote 
of Eq. (4) or Fig. 4. 

The wave vector dependent prefactors, like f q , h q and r q , are strongly coupled to the 
structural input as specified by our model and discussed in the previous chapter. It cannot 
be expected that the simple model reproduces these amplitudes quantitatively. However, 
the variations of the a-amplitude f q c and relaxation time r q c in phase with S q c and of the 
critical amplitude h c q c out of phase with it, are expected to be general findings applicable 
to molten salts. The corresponding variations in the mass-density quantities shown in Fig. 
1 have been found in MCT calculations for different simple liquids [2,6]. They have been 
compared to dynamic light scattering spectra in colloidal suspensions; general agreement 
with errors of the order of 15% was observed [9]. Neutron scattering from CKN measures 
a combination of charge- and mass-density fluctuations determined by the different atomic 
neutron scattering cross sections. Our findings of peaks in the a-amplitudes and in the 
a-times at the wave vectors, q c and q n , characterizing charge and density fluctuations in 
the equilibrium structure, are in qualitative agreement with the reported measurements [7]. 
The results show that charge ordering and local packing are the underlying mechanism. 

Let us restate, that the increase of the time scales at peaks in the static structure 
factors is not a simple example of de Gennes narrowing and cannot be explained by short 
time sum rules but is a consequence of the MCT equations (1) [2]. Eqs. (1) do not include 
the short time dynamics and consequently violate the short time sum rules. The shown 



variations result from the specific magnitudes of the coupling of different modes in the 
memory functions; the couplings are determined by S q . 

The MCT equations (1) correctly describe non-propagating and non-diffusive charge 
fluctuations in the long- wavelength limit (3). This is a consequence of screening in ionic 
melts which requires finite restoring forces for charge fluctuations even on long wavelengths. 
The generalized conductivity changes from a low frequency constant, a oc 1/r, to a power 
law behavior at large lot. 

In a cursory search liquid glass transitions were located for different parameters in this 
model. The parameter of table I lead to a rather strong peak in the charge structure factor, 
Sg C . The origin of this is the large charge asymmetry z 2 /z2 =-3 which entails a corresponding 
concentration ratio due to the requirement of global neutrality. The overestimated charge 
oscillations lead to two special features in our results. First, the maxima in the g-dependent 
amplitudes of the charge fluctuations are quite pronounced. This prevents any quantitative 
comparisons with the neutron or light scattering data of [7,8]. This failure emphasizes that 
quantitative comparisons between MCT calculations for simple liquids and experimental 
data require appropriate microscopic models if non-universal features of the theory are 
tested; see [6,9]. Second, comparing Fig. 1 to the corresponding results for one-component 
liquids, shows that the frozen-in mass-density structure of the formed glass is similar but 
somewhat distorted. Inspection of other parameter values in this model reveals that another 
glassy structure becomes stable if the charge asymmetry is increased some more. The small 
stability of the studied glass to another amorphous structure is the origin of the large 
exponent parameter. This second glass differs in the frozen-in mass-density but not in the 
charge-density structure. A striking consequence of the proximity of the two glassy states 
is the splitting of the a-relaxation in F™ n {t) for some wave vectors into two processes. This 
effect is indicated in Fig. 3, where the inverse a-peak position frequency jumps from one 
process to the other at qa ^2.5. This phenomenon has been discussed in schematic MCT 
models [10] and will be studied further in this microscopic model. 
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